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'—{ ' Abstract. Measured from high-resolution spectra, we present Ca II K line chromospheric fluxes in 35 G dwarf 

stars of 5 open clusters to determine the age-activity-rotation relationship from the young Hyades and Praesepe 
Q>^ • (0.6 Gyr) to the old M 67 (~ 4.5 Gyr) through the two intermediate age clusters IC 4651 and NGC 3680 (~ 1.7 

^SJ ' Gyr). The full amplitude of the activity index within a cluster is slightly above 60 % for every cluster but one, 

NGC 3680, in which only two stars were observed. As a comparison, the same Solar Call index varies by ~ 40% 
during a solar cycle. Four of our clusters (Hyades and Praesepe, IC 4651 and NGC 3680) are couple of twins 
as far as age is concerned; the Hyades have the same chromospheric-activity level than Praesepe, at odds with 
early claims based on X-ray observations. Both stars in NGC 3680 are indistinguishable, as far as chromospheric 
activity is concerned, from those in the coeval IC 4651. This is a nice validation of the existence of an age-activity 
relationship. On the other hand, the two intermediate age clusters have the same activity level of the much older 
M 67 and of the Sun. Our data therefore shows that a dramatic decrease in chromospheric activity takes place 
^Z ' among solar stars between the Hyades and the IC4651 age, in about 1 Gyr. Afterwards, activity remains virtually 

f~^ , constant for more than 3 Gyr. We have also measured v sin i for all of our stars and the average rotational velocity 

shows the same trend as the chromospheric-activity index. We briefly investigate the impact of this result on the 
age determinations of fleld G dwarfs in the solar neighborhood; the two main conclusions are that a consistent 
group of 'young' stars (i.e. as active as Hyades stars) is present, and that it is virtually impossible to give accurate 
Q ' chromospheric ages for stars older than ~ 2 Gyr. The observed abrupt decline in activity explains very well the 

_M I Vaughan-Preston gap. 
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j; 1. Introduction among; which w e mention the Parker's review on the solar 

C^ magnetic fields l|Parkeilll970|) . 

■ - ' Chromospheric activity is powered by rotation in the pres- 
ence of a deep convective envelope. Because of rotation The present study aims at improving this picture 
braking by a magnetised stellar wind, rotational velocity by accurately ascertaining the form of the activity-age 
decreases with age, as well as activity, unless the angu- law on the basis of a new dataset. The decaying law of 
lar momentum is sustained by tidal interaction, as in the chromospheric activity is a precious tool in determining 
case of short-period binaries. This simple picture is at the ages of field main sequence stars and chromospheric ages 
basis of the so called age-activity-rotation relationship in have been used to undertake detailed studies of the star 
solar-type stars and forms the paradigm of stellar activ- formation rate and chemica l enrichment of the Galaxy 
ity in the last 30 years. The hterature on the argument is iJRocha Pinto et al.||2000lj|). 
enormous and the reader may refer to a number of works. The flux emission in the core of the K line of Ca II 

is the most effective tracer to probe the chromospheric 

Send offprint requests to: G. Pace, e-mail: gpace@eso.org , . ., . , , , / rr, 7 71 li ctaI 

^ ^, ■" . „ , , ^r,^ ,,-r^ ^ , activity m solar-type stars (see e.g. Iljlanco et ai.l 119741 

Observations collected at the ESO VLT. Some data pre- ir: — "T^ — :: 1|, ''ij, ^ i i • T" ■ — Tr 

^ , , . u i • J .1 ii 117 Aj T.' 1 r^u M ILmskv fc AvresI 119781 1. Indeed, we intend to use m this 

sented herein were obtained at the W. M. KecJi Observatory, "^ — «v 

which is operated as a scientific partnership among the ^^^^y the equivalent width of the line emission core, cor- 

California Institute of Technology, the University of California, ^ected for the underlying stellar photospheric contribution 

and the National Aeronautics and Space Administration. The and transformed into flux at the stellar surf ace (f^). as in- 

Observatory was made possible by the generous flnancial sup- dex of the level of chromospheric activity l|Pasauini et alJ 

port of the W. M. Keck Foundation. il988i^) . 
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Among similar indices the most used and widely 
known is the Mount Wilson R'fjK^ which is the flux in 
the Mount Wilson photometer passband normaliz ed to 
the stellar eff'ective temperature l|Noves et al.lll984l) . The 
apex indicates that the quantity is corrected for the the 
photospheric-flux contribution. This convention will be 
adopted throughout this paper. While the Mount Wilson's 
index is a contrast index, our F'j^ is an absolute flux. If 
we indicate the chromospheric flux of the H and K line 
cores at the stellar surface (erg cm~^sec~^) as F^^, and 

the bolometric one as Fbol 



then R'uj^ is the ratio ^"'^ . 



J ll99fil:IWilsonlll963t 
relatively little work 



While chromospheric activity has been widely stud- 
icd in field stars (sec e.g. Henrv et 
LPasquini et al. 1990: Pasquini 199S 
has been done on clusters, which are the perfect calibra- 
tors of the age-activity-rotation relationship: they are in 
fact homogeneous samples of stars with well determined 
ages. 

For the reas ons explained in 

Ivan den Bergh fc McClurj (|l980(), old open clusters 
are rare an d tend to be distant. This is why the original 
ISkumanichI l|1972f) law was based only on 1 point older 
than the Hyades, which was the Sun. Only with the ad- 
vent of 4-m telescopes and high-sensitivity instruments 
has it been possible to acquire Call spectra for subgiants 
in M 67 tiDupree et al. 1999). wh ich is as distant as 850 
pc (see e.g. ICarraro et al.l 1199(11) . Nowadays 8-m class 
telescopes are sufficiently sensitive to allow the study of 
Ca ll emission even in the dwarfs of distant clusters. 

iBarrv et al.l l|1987J) made the most consistent attempt 
to cover the lack of old-intermediate age clusters data. 
They built a 2-A resolution system and used it to observe 
several intermediate and old clusters. This very interesting 
attempt, on the other hand, needs to be looked with some 
caution, since, for instance, the typical G star emission 
core is much larger than 2 A. Therefore a 2-A resolution 
is sensitive to other factors, such as zero point correc- 
tion s and photospheric profi le in addition to activity (see 
also ISoderblom et al.lll99ll for a discussion on Barry et 
al.'s data). A narrower, pure chromospheric-activity in- 
dex based on high-resolution spectroscopy, such as the 
quoted F'j^, is of course much more reliable. This is the 
main ground of the present work. 



2. Observations and data sample 

Our study is based on spectra of 3 5 stars in 5 open clusters 
and on the solar data from Whi te fc Livingston (1981) . 
White and Livingston made a very detailed study of the 
Call solar K line covering the cycle 21 from the first mini- 
mum to the maximum. Because of differences in data anal- 
ysis and reduction, their 1-A indices have to be multiplied 
by 1.11 to be homogeneous with ours. 

The Praesepe, NGC 3680, IC 4651 and M 67 observa- 
tions were obtained with the UVES spectrograph at the 
VLT Kueyen telescope. The blue slit was set at 0.7 arc sec- 
onds giving a resolving power of R ~ 60000 ( Dekker et all 



|200(]|) . For the faintest stars, up to 3 consecutive 90- 
minutes exposures were co added, in order to reach enough 
signal at the bottom of the deep Call K line which is typ- 
ically less than 10% of the continuum in non active stars. 
The signal to noise ratio per pixel at the bottom of the K 
line ranges from about 10 for the worst spectra in M 67 
to about 30 for most of the spectra in Praesepe. 

The red slit was set to 0.4 arc seconds, to obtain a 
resolution of R=100000 

The data were reduced using the UVES pipeline 
iJBallester et al.l200(]|) . radial velocity corrected and finally 
co-added. 

The Hyades spectra were obtained within the program 
of search for extra-s olar planets among the solar-type 
stars of this cluster l|Cochran et alJ l2002t IPaulson et alJ 
120021) with HIRES at the Keck Observatory. They have a 
resolution similar to the blue UVES ones (R=60000), and 
a signal to noise ratio at the bottom of the K line ranging 
from about 20 to about 30. The original Hyades' sample 
comprises a large range of spectral types, and was selected 
to include only stars with t;sini lower than 15 km/sec. 
We have selected for the present study only the dwarfs 
with B— V colour s simil ar to the solar one. According to 
IStrassmeier et al.l f2000") (see the g panel in their Figure 
4), late-type stars which rotate faster than 15 km/sec are 
rare exceptions (unless we ar e dealing with clusters much 
younger than the Hyades, see lBouvieilll997J) . therefore we 
do not expect to have a strong bias towards slow rotators. 

All the Call data have been norm alized to the 3950.5 
A pseudocontinuum ()Catalanolll978l) . 

The clusters span a wide range of ages, from the young 
Hyades and Praesepe (0.6 Gyr) to the old M 67 (4.8 Gyr). 
Together with the Sun, they constitute three pairs of ob- 
jects at three different ages: Hyades and Praesepe, IC 4651 
and NGC 3680 (1.6 and 1.8 Gyr), and M 67 and the Sun 
(4.8 and 4.48 Gyr). This allows to have a first check on 
the consistency of the age-activity-rotation relationship: 
if each pair has the same chromospheric-activity level, we 
will have a strong confirmation that chromospheric activ- 
ity is an age driven parameter. 

In addition, the comparison between Hyades and 
Praesepe is extremely in teresting, after the finding by 
iRandich fc Schmitti l|1995f) that the solar stars belonging 
to these clusters have different X-ray (coronal) emission, 
with the Hyades showing sub stantially higher X-ray lu- 
mino sities than Praesepe f see iBarrado v Navascues et alJ 
Il997l and reference therein). 

T he ages of M 67, N GC 3680 and IC 4651 are taken 
from ICarraro fc Chiosil l|1994'l . Hyades and Praesepe, 
with some differences between various authors, are often 
thought to be coeval with ages of about 0.6 G yr (se e e.g. 
IPerrvman eTLx. 1998: Mermilliod 1981: AlleJfl973() . For 
the solar age we have used Guenther's work (1989), based 
on the age of the oldest meteorites. 

All the stars observed with the VLT were chosen to be 
on the main sequence, high-probability members of the 
clusters and not known to be binaries at the moment of 
the observation. The choice was done using as references 



G. Pace and L. Pasquini: The Age-Activity-Rotation Relationship. 



iNordstrom eHIl t99(i] for NGC 3680, iMeibom et all 
ll2002l) for IC4651 and s everal sources for M 67, includ- 
ing |LathaineL§iJ l|l992(l for the binary determination in 
this cluster. Our stars have colours comprised within 0.51 
<(B— V)o < 0.73, closely encompassing the solar colour 
which, with some variations between different authors, is 
evalu ated around B— V=0.63 (see e.g. IStraizis fc Valiugal 
I1994I) . We also note that the paucity of stars observed 
in NGC 3680 is i mpose d by the fact that, as shown by 
INordstrom et al.1 <|l997|^ . this cluster has very few single 
G stars members left, most of its low mass members be- 
ing dispersed in its lifetime. 

As far as Praesepe is concerned, out of the 7 
stars observed we c hose 3 stars detected in X-ray by 
iRandich fc Schmitti l|l995ll and 4 undetected, in order to 
check whether the chromospheric activity was different be- 
tween the two groups, as expected if chromospheric and 
coronal activity were strongly correlated. 

All the observed stars are stimmari zed in Table ^ 
Star names are taken f r om | Eggenl 1119691) for NG C 3680, 
lAnthonv Twarog et all l|l988l) andlEggenl lll97lh for IC 
4651. ISande rs' 1*1977^ for M 67, 'Klein WassinkI l| l927ll) for 
Praesepe and van BuerenI (|l952j for the Hyades stars. 
The photometr y is taken from the refe rence papers for 
IC 4651, from Nordstrom ct al.' (l99t) for NGC 3680, 
from Hipparcos Catalog ue (ESA 1997'! for Hyades, from 
Ijones fc Stauffed l|199lj) for Praesep e. with the excep- 
tion o f KW 368 which is taken fro m Jone s fc CudworthI 
l|l983(l . and from lMontgomerv et al. U993) for M 67. For 
Hyades no corre ction at all for inters tellar (IS) extinc- 
tion is applied JTavlor fc .Tonal l2002|) . and NGC 3680 
colour s are given already dereddened bv INordstrom et al] 
|l993). B-V cofours of Praespe, M 67 and IC 4651 had 
to be corrected for colour exce ss. We used E(B— V )=0.05 
for M 67 (also taken fror n iMontgomerv et al.l Il993i). 
_E(B- V)=0.086 for IC 4651 iJAnthonv Twarog fc Twarog 
19871) and E(B-V)=0.02 for Praesepe iJTavlor fc .Tone: " 



20021) . The uncertainties in reddening are not at all critical 



for the analysis performed. 

The photometric data retrieved from the literature are 
shown in the second and third columns of Table ^^V mag- 
nitude and B— V colour index, respectively). The differ- 
ent chromosphcric-activity indices are reported in Table 
^ from the fourth to the sixth column and explained in 
the caption. 

The quantity shown in the seventh Column is the Full 
Width Half Maximum (FWHM) of the cross correlation 
profile of the spectra. Th e cross corr e lation profiles were 
obtained as described in iMelo et alJ l)200lj) : we used as 
template a box-shaped mask, suitable for solar-type stars, 
which was kindly provided by Claudio Melo. We processed 
the red part of the red arm UVES-spectra, all rebinned 
at the same step (0.0174087 A). The la uncertainties of 
the FWHMs are of about 1%. 

The rotational velocities determined by us are given in 
the last column . They are obtained from the FWHMs via 
the calibrations of Section rOl 



It is worth noticing that what really matters for the 
present purpose is the evolution of the angular momen- 
tum of the stars, we are therefore mostly interested in the 
differential vsini values. 



3. Analysis 

3.1. Call fluxes 

The analysis is based on the measurement of flux in the 
emission core of the Call K line. This line is the most sen- 
sitive indicator of chromospheric activity, and it has been 
widely studied in a number of dif ferent environments, in- 
cluding the Sun a n d solar stars 
1 1 I r 



Linsk^^^al. 1979; Pasqui ni et al 



S Linskv fc Avretd 
119881: IPasauinil 



197C 



1992 



White fc Livingston.l98D . 



In general, the most sensitive area is that comprised 
between the line Kl minima, but those are hardly seen in 
low activity stars or in intermediate S/N ratio spectra. 

A very useful index of the level of chromospheric ac- 
tivity is the area of the 1-A wide band centered in the line 
core (hereafter referred to as 1-A index). This 1-A band 
corresponds roughly to the separation between the Ki 
minima for an active star (as those in Hyades or Praesepe), 
and it is about 30 % larger than th e Ki minima separa- 
tion f or a star as act ive as the Sun l|White fc LivingstonI 
Il98ll: IPasauinilll992l ). The 1-A index was measured in 
the spectra either uncorrected or corrected for IS absorp- 
tion. The results are shown in the fourth and fifth col- 
umn of Table H respectively. In addition, we recall that 
the F'j^ is a measurement of a contrast between the chro- 
mospheric emission and the underlying photospheric flux. 
Several calibrations can be used to eliminate this offse t 
and to quantify different band widths llLinskv et alJll979l) ; 
this is widely discussed (e.g. in lPasouini et alJll988ll . and 
will not be repeated here. 

In addition, although our clusters are among the clos- 
est and less reddened, the Ca II lines could be contami- 
nated by the IS Call, a possibility that we have taken into 
account as explained later in this Section. 

We summarize the different steps required to get the 
final fluxes from the observed spectra. 

1. Spectra normalization and rcctiflcation 

2. Subtraction of the IS line 

3. Index measurement 

4. Subtraction of the photospheric contribution 

5. Calibration in Fluxes 

The observed spectra have been n ormalized to the 
3950.5 A ps eudocontinuuni . as done in ICatalanol 1)19781) 
and then in IPasauini et alJ Ijl988') . Since the spectra are 
the result of a merging of different echelle orders, we must 
check that merging is properly made. We did it in two 
ways: by observing early-type stars (i.e. stars with no Ca 
II K stellar line), and by comparing the target spectra with 
a synthetic solar spectrum (courtesy of P. Bonifacio). This 
comparison of observed with photospheric-line profiles has 
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been largely used by us, not only to verify that the rectifi- 
cation of the spectra was properly done, but also to check 
that other effects, such as background light, were correctly 
taken into account. The K line is very suitable for this 
purpose, because it is very deep, and therefore sensitive 
to small mistakes in setting the zero flux level, fn Figure 
Q]we show the average spectra of all the cluster stars and 
the solar one, all plotted over the synthetic solar Call K 
stellar line. The solar spectrum is taken from the UVES 
archive ^ . We checked that the observed photospheric-line 
profiles successfully overlap after normalization. Of course 
this is only possible because all of the stars have similar 
colour to that of the Sun. We have also checked the con- 
sistency of each observed K line profile with the synthetic 
solar K line profile. The spectra of the single stars are 
shown in Figures |21 and |3| 

For all clusters only one component of the IS Ca II 
K line falls in the stellar-spectra region where the 1-A 
index is measured. In order to evaluate its contribution, 
we have observed for each of the 3 reddened clusters a 
reference star with an easily recognisable IS component. 
These stars show no narrow stellar line inside the Ca II 
absorption feature either because they are fast rotators (as 
in the case of Sanders 1082 in M 67, an Algool-type binary, 
and Eggen 70 in IC 4651) or because they are hot enough 
(as in the case of Eggen 10 in NGC 3680) to present a 
pronounced broad deep core. In Figure 01 upper panels, 
the spectra of the quoted reference stars are shown. For 
each of them, the spectrum without the IS contribution is 
plotted over the original one. The latter has been divided 
by the former, and the result is the IS contribution, which 
is shown in the bottom panel of Figure^ 

Eggen 10 is quite hot, but this is not the case for 
Sanders 1082 and Eggen 70. If they were slightly active, 
their emission would partially fill in the absorption IS 
line. In this case, the IS negative contribution would be 
underestimated by the removal procedure. We are, how- 
ever, confident that our estimate is quite robust, because 
Sanders 1082 and Eggen 70 are rather fast rotators, there- 
fore if their chromospheric emission were relevant, they 
would be as broad as the photosphcric lines and the emis- 
sion wings would be visible outside the IS contribution. 

Each K line spectrum of our target stars has been di- 
vided by the appropriate IS spectrum, and the 1-A indices 
measured again. The results are given in the fourth col- 
umn of Tabled In FigureElwe plot the F^ vs. B— V before 
and after the IS-line contribution has been subtracted (re- 
spectively left and right panels). Clearly the effect of this 
subtraction, although noticeable, does not affect dramat- 
ically the relative emission between the different clusters, 
which is not surprising, given the low reddening of the 
clusters. 

Chromospheric flux es decrease for cooler stars (see e.g. 
ISoderblom et al.l 119911) . In Figure El such a trend is not 
present, with the possible excepti on of Praesepe stars. By 
comparing our plot with that of ISoderblom et al.l l)199ll 

^ http://www.eso.org/observing/dfo/quaIity/UVES/pipeIine/ 
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Fig. 1. Average clusters Call K line spectra. Overimposed 
is the synthetic photospheric-solar model. In the bottom 
panel we also plotted the real solar spectrum, also overim- 
posed on its model. The effect of the back heating of the 
chromosphere is evident also outside the core. 



see Figure 6 therein), we do not expect our results to be 
significantly affected by the different colour distribution 
of the targets in the different clusters. 

We have subtracted to the 1-A indices the pho- 
tosphcric contribution, calculated on the synthetic 
photospheric-solar model. The quantity subtracted 
amounts to 0.049 A. Note that, even though this sub- 
traction is the same for all of the stars, once the 1-A 
indices have been transformed into fluxes in the way ex- 
plained later in this Section, the photosphcric contribution 
becomes colour dependent, because the transformation 
into flux is colour dependent. The photosphcric contri- 
bution we subtracted does not differ more than 20% from 
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Fig. 2. HIRES spectra of the Hyades stars of our sample. 



other appro ximations, e.g . the o ne used by iNoves et alJ 
l)l984 - and iDuncan et all l)l984l) . The choice of another 



approximation does not significantly change the result, 
even though it might affect the spread within each cluster. 
We have derived the chromospheric fluxes at the stellar 
surface transforming the IS corrected, photospheric sub- 
tracted equivalent widths into flux densities at the stellar 
surface (erg -sec"^- cm^^). This was made in three steps. 
We first computed the V— R colours from the B— V ones. 
We used the calibration in lPasauinil (|l985J): 

V -R = 0.618 • {B - Vf - 0.183 ■{B-V) + 0.375 (1) 

We have then used the calibration in lPasauini et alJ l|1988j) 
between the V— R colour and the flux at the stellar surface 
of the pseudo-continuum radiation at 3950.5 A: 



logF. 



3950.5 



8.459 - 2.833 ■ {V - R) 



(2) 



We have finally multiplied the equivalent widths for 
-F3950 5 , obtaining the fluxes in the sixth Column of Table 

m 

We need to estimate the absolute and relative accu- 
racy of the indices and of the fluxes. Actually, the relative 
uncertainties are here the most relevant, since we will ulti- 
mately compare the ratio of the fluxes between the differ- 
ent clusters. The absolute value, although interesting, is 
just a normalised value for the purpose of the age-activity 
law. 

The Signal to Noise per pixel is not the limiting fac- 
tor of the accuracy, even though it is only of '-^ 10 in the 
spectra of the faintest stars, because we integrated on 1 A, 
namely ~ 48 pixels and ^ 37 for HIRES spectra) . For the 
intermediate age clusters and M 67, the dominant sources 
of error in the measurements could come from an incor- 
rect treatment of the inter-order (and scattered) light or 
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Fig. 3. UVES spectra of the stars of our sample. For the stars in IC 4651, NGC 3680 and M 67, in each panel, the 
spectrum non corrected for IS line is plotted over the corrected one. 
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Fig. 4. Template spectra used to detect the IS-line feature. The original spectra (one for each of the 3 distant clusters) 
are over plotted on the supposed non contaminated profile on the 3 upper panels. In the lower panels the resulting IS 
features are shown. 
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Fig. 5. The trend of F^ with the colour. In the right panel the chromospheric activity is not corrected for the IS 
absorption, in the left one chromospheric-activity indices regarding IC 4651, NGC 3680 and M 67 stars are measured 
on IS absorption corrected spectra. There is no remarkable difference between the two diagrams. A decreasing trend 
with colour for chromospheric-activity index can only be seen for Praesepe stars, since the range of colours involved 



IS narrow. 



from rectification problems. As we have mentioned before, 
we think we are rather free from these problems, because 
of the satisfactory match between the solar synthetic and 
stellar observed profiles of the K line wings (i.e. the region 
between ~ 3920 and ~ 3950 A, outside the line core). We 
have first divided the photospheric spectra of all the stars 
of each cluster for their average spectrum, finding peak 
to peak differences of less than 7% in their profile ratios. 
We have then repeated the same comparison with the solar 
synthetic spectrum; and the differences between the stellar 
spectra and the solar model are greater. They arrive up to 
20 % for Hyades and Preasepe stars, while they are of up 
to about 10% for the other clusters. Hyades and Praesepe 
stars show exactly the same behaviour, despite the spec- 



tra were obtained with two different spectrographs. The 
stellar spectra have always higher profiles than the solar 
model and that is true even for the Sun itself (cfr. Figure 
P). This can be due to a low level of residual scattered 
light in the spectrograph, but the main reason is that the 
solar model we used does not include any chromospheric 
effect, which will tend to heat the higher layers of the pho- 
tosphere and therefore to raise the photospheric profiles 
as well. This effect is well known at least in the Sun, where 
the differences between the quiet and plage re gions arise 
also in the photospheric wings of the K line I La Bontd 
Il98a. see Figure 1 therein). 

Since the IS correction could be rather tricky, we have 
tried several different fits of the IS-line wings finding sub- 
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stantially only minor variations and we estimate the error 
from a wrong computation of the IS line as negligible. It is 
important to see also that this error is the same in both, 
absolute and relative terms and that the whole IS correc- 
tion accounts for at most 20% of the chromospheric flux 
measured. 

The errors in F^ come from three main sources: 

- the uncertainty in the equivalent width measurement, 
which on its turn has, as mentioned, many other 
causes; 

- uncertainty in the V-R vs. B-V calibration; 

- the uncertainty in the B-V colours. 

The absolute flux determination has an error of 30%, 
due to the uncertainties in the logi^agso.s vs. {V — R) cal- 
ibration. However this affects mainly the zero point of our 
F'j^ measurements, and only to a much smaller extent the 
relative comparison. 

The two other sources of error give a contribution to 



the relative uncertainty 



F' 



that can be calculated on 



the basis of Equations ^ and |5| They are, respectively: 

AEW 



and 



EW 



(8.05 ■(B-V)- 1.19) ■ A(B - V) 



The typical relative maximum error in the equivalent 
width measurement of the single spectrum is within 10%. 
For an error in the colour A{B — V) = 0.04, which we 
adopt as value of the maximum error in the colour, the 
corresponding relative error in the flux is, depending on 
B-V, between the 12 and the 18%. If we add quadratically 
the two contributions, we flnally obtain a maximum error 
in the single star measurement ranging from the 15 to the 
20%, if only relative values are concerned. 

3.2. Determination of the vsini 

As stated in Section [21 we have measured the FWHM 
of the cross correlation proflle for all the spectra of our 
sample, and the data are shown in the seventh column of 
Table m 

Lin e broadening, as extensively discussed e.g. in lGravl 
Ijl992l) . is caused by several effects - such as thermal mo- 
tions, micro- and macro-turbulence. Stark effect, Zeeman 
effect, natural broadening, instrumental response and ro- 
tational velocity. We make two assumptions. The first is 
that the sum of these effects produces a gaussian shape, 
whi ch is a reasonable as sumption for solar-type stars (see 
e.g. lPasauini et al.ll988l) . The second is that the difference 
in broadening among the stars of our sample is due only to 
rotation, or that the quoted effects contribute equally to 
the line broadening for all the stars, which is well justified 
given the similarity in the spectral type of the targets. We 
can now write: 



where uccp is the sigma the cross correlation profile, 
provided that it is a gaussian. For the reasons explained 
above, CTq is expected not to change significantly from star 
to star, within our sample. Since <Jccp and cr^sini are pro- 
portional respectively to the FWHM of the cross correla- 
tion profile and to vsini, we expect vsini to be: 



usmi 



i = A- VfWHM^ - B2 



(3) 



where A and B, like ao, are constants. Using vsini mea- 
surements from the literature we can calibrate A and B in 
the previous relationship, and then calculate vsini for all 
the stars. We have done this using the data of Table [3 in 
particular the v sin i measurements in the third Column 
and the FWHMs in the fourth Column. 

Most of the spectra of the calibrators have been ob- 
taine d in a parallel project on IC 4651 i Pasauini et alJ 
|2003|) . The spectrum of the Sun is from the UVES archive, 
and is the same plotted in the bottom panel of Figure ^ 
For Eggen 60 we used the FWHM of Table ^ since it is 
one of the targets. All the spectra used for the calibra- 
tion are taken with the same instrument (UVES) at the 
same resolution (R=100000) of our observations and the 
same spectral range is used to get the cross correlation 
profiles. The vsini are taken from|Mgi^22L£ii^yi 
for the stars in IC 4651 and from lNordstrom et al.l 
for Eggen 60. For the Sun we adopted 2 Km/sec. 

We decided not to use for the calibration Eggen 70 in 
NGC 3680, even though for this star vsini is available. 
The reason is that the 1 a error given for v sin i is three 
times larger than iisini itself, that is 1 km/sec. 

The calibration is a x^ fit to the data of the Equation|31 
The root mean square of the residuals of the data around 
the fitting curve is 0.26 km/sec, no star deviates more than 
10% of its vsini, and 5 out of 8 calibration stars deviate 
less than 2%. 

The result of the calibration is: 




; sini = 0.72 • {FWHM^ - 9.05^ 



,1/2 



(4) 



f^CCP 



— \J dv 



sini^ + CTo^ 



The curve is plotted in Figure |B| 

We have used it to compute the values given in the 
eighth column of Table ^ We remind the eventual users 
that we used the spectral range 588-680 nm, a spectral 
resolution of R= 100000, and a rebinning step of 0.0174 A 

We note that the small residuals from our x^~fit pos- 
sibly overestimates its real accuracy. The reason is that 
our calibration is a secondary calibration, because also 
the calibrating measurements are based on the FWHM 
of the cross correlation profiles. The errors may be quite 
larger in particular below v sini ^ 3 Km/sec, where other 
causes of line broadening become more important. We un- 
derline that we are mainly interested in differential values 
of vsini, that are as reliable as the root mean square of 
our calibration indicates. 

Only for Hyades stars we could not use this calibration, 
since their spectra were taken with another spectrograph 
(HIRES) at a different resolution. Anyway, all of the v sin i 
for our Hyades targets are measured in IPaulson et alJ 
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I|2003j) . In this work, the rotational velocities are deter- 
mined by comparing five synthetic Fe I lines with the ob- 
served ones, therefore temperature is taken into account. 
In Figure [7| we plotted the diagram wsini vs. FWHM. As 
for the calibration of UVES spectra, we computed a fit of 
an equation of the type of Equation |31 to the data points. 
The result is that only one star deviates from the fit more 
than 5 % of its v sini, and the median deviation is of the 
2%. This proves that, for targets comprised in the short 
B— V colour range 0.51-0.73, the FWHM is virtually com- 
pletely determined by usini. The presence of an outlier, 
namely van Bueren 88, confirms that more care has to 
be exercised in dealing with slow rotators. We will not 
use the outlier value of v sin i in the analysis of the angu- 
lar momentum evolution. On the other hand, the overall 
quality of the data is evident when comparing Paulson et 
al.'s V sin i values with equatorial velocities from rotational 
periods. 

iRadick et alJ Ijl987r provide precise period determi- 
nation of 5 of our targets. They used narrow band by 
photometry to obtain light curves of 24 Hyades members. 
From the light curves it is possible to obtain the rotational 
period of the stars because rotational modulation of pho- 
tospheric features (spots) induces photometric variations. 
We transformed the periods into rotational velocities, and 
found the data in Table 01 We had to estimate the stel- 
lar radii, which were computed using Tg// and computed 
bolometric luminosity. We transformed th e B— V colour s 
into Teff and bolometric corrections using iFloweJ l|l996l) , 
we then used Hipparcos parallaxes l|ESAlll997r to trans- 
form the apparent magnitudes into luminosities. 



van Bueren Period inferred v 

IQ. [days] [km/scc] 



31 


5.45 


11.1 


52 


8.05 


7.4 


65 


5.87 


10.5 


73 


7.41 


7.9 


97 


8.55 


6.2 



Table 3. Rotational periods of Hyades from lRadick et alJ 

(jl987) and rotational velocities from them inferred. The 
data are collected to compare equatorial with projected 
rotational velocities. 



SUN «>Eggen 




10 20 30 

FWHM [pixels] 



STAR 



Cluster 



V sin i ]km/scc] FWHM [pixels] 



AT 1228 


IC 4651 


15.7 


Eggen 15 


IC 4651 


10.0 


Eggen 34 


IC 4651 


25.2 


Eggen 45 


IC 4651 


4.2 


Eggen 79 


IC 4651 


21.8 


Eggen 99 


IC 4651 


28.1 


Eggen 60 


NGC 3680 


1.9 


Sun 


- 


2.0 



22.24 
16.65 
37.18 
10.27 
31.42 
40.45 
9.89 
9.41 



Table 2. Literature data used for the calibration in Figure 
El The range of B— V colours is 0.514-0.641, closely encom- 
passing the solar one. 



We overplotted the rotational velocities of Table (31 on 
the diagram in Figure d They are equatorial velocities, 
i.e. upper limits of the projected rotational velocity vsini, 
which is the quantity really related to the FWHM of the 
cross correlation profile. On the other hand, the fact that 
rotational modulation is observed indicates that the incli- 
nation angle of stars cannot be too small. We therefore 
expect rotational velocities from period determinations to 
be right where they are: slightly above the projected ro- 
tational velocities. 



Fig. 6. Calibration of vsini with the FWHM of the cross 
correlation profiles of UVES spectra. 

4. Discussion 

After showing in the previous sections that we were able 
to derive chromospheric fluxes and rotational velocities 
with great accuracy for our cluster stars, we can now pro- 
ceed to the discussion of the results. We first investigate 
the spread in the level of chromospheric activity in each 
cluster, which will tell to which extent a single value can 
be used for each cluster. By comparing the couple of twin 
points we can then evaluate if chromospheric activity is 
a function of age for solar stars. Only after that we will 
be able to derive an age-activity relationship. Finally, the 
same is applied to rotational velocities. 

4.1. Variability within each cluster 

Accor ding to the detailed study of IWhite fc Livingston! 
lll98l|) . the level of solar chromospheric activity varies 
of ~ ±10% (from peak to peak) during the 11-years 
solar cycle. If we consider the data corrected for pho- 
tospheric contribution, then the variation becomes of 
^ ±20%. Analog cycles take also place in solar-type 
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[pixels] 



Fig. 7. wsini vs. FWHM diagram for Hyades stars, 
and the fit to the data points: vsini = 0.82 • 
[FWHIvP - 10.75^) . Also equatorial velocities from 
periods of rotation are plotted. The latter are indicated 
as square symbols with the downward arrow (they are in 
fact upper lim i ts for tisinz), while wsini data points from 
IPaulson et al.1 1)2003() are indicated as triangles. 



stars (seee.g.lHempelmann et al.ll996HRadick et al.ll998t 
iBaliunas et alJll99^ . This poses some caution in deriving 
chromospheric ages from a single spectrum, and causes 
scatter around the mean chromospheric-activity levels ob- 
ser yed in each cluster, i.e. at a given age. 

iGiampapa et al.l 1)2000(1 have shown that indeed in the 
old M 67 considerable variations are observed if a sample 
large enough is acquired: they find that approximately 
30% of the solar-type stars in M 67 exhibit levels of ac- 
tivity that are outside the envelope given by the solar- 
activity cycle. 



cluster 



spread stars 



PRAESEPE 


59.3% 


6 


HYADES 


68.3% 


15 


IC 4651 


60.2% 


5 


NGC 3680 


31.2% 


2 


M 67 


62.4% 


7 



SUN 



43.7% 



Table 4. Peak to peak spread of the F'^^ fluxes within 
each cluster, expressed in percentage of the related mean 
value. The spread of the solar data is the variation during 
the solar cycle 21. 



As shown in Table |2 all clusters but NGC 3680 
(in which only two stars are observed) have a peak 



to peak spread shghtly exceeding ±30% of the mean 
chromospheric-activity level. For young clusters such as 
the Hyades and Praesepe we expect that part of this vari- 
ation is due to the superposition of long and short-term 
variability, the last induced by rotation. Another huge 
cause of chromospheric variability in young stars is flar- 
ing activity. For M 67, however, for which we expect the 
members to behave much like the Sun, the scatter is 1.4 
times larger than the variations during a solar cycle. This 
means that an additional spread of ± 24 % has to be ex- 
plained. For M 67, whose mean B-V colour is 0.62, the 
measurement maximum errors are of about the ± 18 %. 
This may indicate that cycles similar to the 11-years solar 
cycle do not fully account for the dispersion of the data, 
with interesting implications for the study of the Earth- 
Sun system. Admittedly, we have observed each star only 
once (three consecutive times the fainter stars), rather we 
assume that the stars in our sample have been observed at 
only one phase of their long-term activity, which is likely 
different from star to star. This implies, for instance, that 
possible star to star differences may cause additional noise 
to our results. We think, however, that our procedure of 
overlapping the photospheric proflles eliminates most of 
the observational uncertainties. Before inferring firm con- 
clusions, more data have to be collected. 



4.2. Hyades vs. Praesepe, IC 4651 vs. NGC 3680, M 
67 vs. Sun 

The Hyades and Praesepe clusters are particularly rele- 
vant because they have b een studied also in the X-ray 
IIRandich fc Schmitd Il995j) . It is at least reasonable to 
reckon that the X-ray emission (i.e. coronal activity) and 
the chromospheric activity are not completely indepen- 
dent, since t hey are both related t o rotation and magnetic 
field fs ee e.g.lBruevich et al.l200l[ and references therein), 
but in lRandich fc Schmitd fl995j the two coeval clusters 
were fou nd to have substantially di fferent levels of X-ray 
emission ( Randich fc Schmittil 199511. This early result has 
been recently revised bv lFranciosini et alJ ()2003ll . who find 
that Praesepe stars have the same level of X-ray activity 
of Hyades if the central field is considered. 

Our data confirm that Praesepe and Hyades stars 
span the same ranges of chromospheric-activity levels and 
V sin i values, even though our Praesepe stars have a mean 
X-ray luminosity significantly lower than the Hyades ones 
(less than a half). 

In Figure IHl the Lx vs. vsini (right panel) and the 
Lx vs. F^ diagram (left panel) are shown fo r both clus- 
ters. Lx measurements are taken from fro m IStern et alJ 
11995) for the Hyades, and [Randich & Sc hmitd 1)1995^ 
for Praesepe mernbers. but KW 326, which is given in 
iFranciosini et alJ l|2003>l 



In both panels the spread is high and no real relation- 
ship is observed. In spite of the aforementioned possible 
discrepancy between X-ray emission and chromospheric 
activity, a significant correlation likely exists. The fact 
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Fig. 8. X-Ray luminosities vs. chromospheric activity and rotation diagrams. In neither of those a trend can be clearly 
seen. 
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Table 5. X-ray flux measurements of Praesepe and 
Hyades stars. All of them but one refer to the ROSAT 
detection in the band 0.1-2.0 keV. KW 326 was not de- 
tected by ROSAT. Its measurement, which refers to the 
same band, has been performed with XMM-A^ewton satel- 
lite. 



that it is not seen in Figure |2| could be due to the ac- 
tivity variability, since the optical and X-ray observations 
are not simultaneous. 



When focusing on the other coeval clusters we notice 
in Table^that the activity and rotation values for IC 4651 
and NGC 3680 are, as in the in the case of Praesepe and 
Hyades, indistinguishable. We have to remind the reader 
that only two stars were observed in NGC 3680, there- 
fore its average chromospheric-activity and rotation lev- 
els cannot be considered accurate. Nevertheless it is very 
significant that both NGC 3680 stars have values similar 
to those of IC 4651. This result is particularly important 
because, as we shall see in the following, these clusters 
show a novel behaviour with respect of what thought so 
far. 

Finally, given the excellent agreement between the 
mean M 67 activity values and the Sun, we can conclude 
that all our coeval points show very consistent average 
values, therefore the level of chromospheric activity is an 
age driven parameter. 

The V sin i for M 67 is slightly larger than the solar one 
{^ 3 km/sec, see Section IT^ : we do not consider this at 
the moment as a significant discrepancy, but it would be 
interesting to verify that M 67 cluster stars rotates indeed 
slightly faster than the Sun. 

4.3. The age-activity-rotation relationship. 

When the stars of about 1 Mq reach the ZAMS, they show 
a high dispersion in rotational velocity, and therefore in 
chromospheric-activity levels, as a result of their differ - 
ent Pre Main Sequence histories (see e.g. lBouvieiill997l) . 
But at 600 Myr, i.e. the youngest age we are dealing with 
in the present study, the magnetic braking mechanism al- 
ready made the rotational velocity uniform. So, for each of 
our clusters, rotational velocity and chromospheric activ- 
ity are expected to be all of the same order of magnitude 
(see Section ^TJ. The evolution of these quantities, as it 
results from our data, is shown respectively in the two 
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panels of Figure |^ It is very important to bear in mind 
that in this Figure the bars are not error bars, but they 
represent the peak to peak spread within the clusters, i.e., 
the difference between the maximum and the minimum 
cluster values, and the points are the mean values. 

As expected, chromospheric activity and rotation show 
the same decaying trend. Surprisingly, however, the inter- 
mediate age clusters IC 4651 and NGC 3680 have the same 
level of activity and rotation as the older M 67 and the 
Sun. We recall that this two data points are based on mea- 
surements of 7 stars, two in NGC 3680 and 5 in IC 4651, 
none of which shows a chromospheric activity above the 
most active star in M 67, that is to say that the activity 
distribution of the 7 intermediate age stars and that of 
the 7 M 67 stars are virtually equal. 

This means that an abrupt decay of chromospheric 
activity occurs between Hyades' age and 1.5 Gyr, after 
which a very slow decline follows, if any. Three power 
law decay laws are plotted over the data points of Figure 
El the two more gentle slopes re present the t~^^'^ and 
^-2/3 iQT^g proposed resp ectively bv lSk^imanichl Ijl972l) and 
ISoderblom et al.l l)l99l|) . scaled to best match our data 
points. 

The intermediate age clusters are clearly not consis- 
tent with any of the laws so far suggested; they have an 
activity level lower more than a factor 2. If the exponent 
of the power law is left as a free parameter, we find that 
the points are best fitted by the third curve plotted on 
the F'l^ vs. age diagram, namely Fk — 0.5 • t~^-^ + 0.5, 
where F^ is in unit of 10^-erg-cm~^-sec~^, and the age is 
expressed in Gyr. Note that neither a different zero point, 
nor different assumptions on the age of the clusters would 
change our conclusion, since activity keeps constant for a 
time interval, the one between the intermediate age clus- 
ters and M 67, that under no acceptable assumptions can 
be shorter than 2 Gyr. Similar results are also claimed by 
iRodrigues et all l|2003^ . 

Figure recalls what observed among field stars: the 
so called Vaughan-Preston gap, i.e. the lack of nearby 
stars with a level of chromos pheric activity interrnediat e 
between Hyades' and the Sun ijVaughan fc Prestonlll98(ll) . 
With our results this gap gets a quite natural explana- 
tion: the decline in activity is so abrupt to be almost a 
discontinuity, and only stars with ages in the range be- 
tween 0.6-1.5 Gyr may fall in it. 

As mentioned above, the previous most systematic at- 
tempt to study chromospheric activi ty in stars be l ongin g 
to open clusters has been made by iBarrv et al.l ijl987l) , 
who obtained results quite different from ours. Some ma- 
jor discrepancy may exist i n the measurement of the ac- 
tivity index. For instance, in' Barrv et a l. '1987), the ratio 
between the chromospheric fluxes in the Hyades and in M 
67 stars or in the Sun, is only ^ 2.2, while for us is larger 
than 4. We find a similar ratio (~ 2.4) if our index includes 
the photospheric contribution which represent a fraction 
of the 1-A index much larger for M 67 than for Hyades 
stars. In Figure|n|thc difference between ours and Barry et 
al.'s data is apparent. Given the large scatter of fluxes ob- 



served in each cluster, we cannot, however, exclude that 
the discrepancy is due to differences in the sample ob- 
served. We tend to consider this possibility unlikely, con- 
sidering the large number of stars observed. Barry et al. 
have observed 13 stars in the Hyades and 16 in M 67 and 
we have observed 15 stars in Hyades and 7 in M 67. These 
authors measure also a chromospheric-activity index for 
the intermediate age cluster NGC 752 only 1.4 times lower 
than the Hyades and 1.5 times bigger than M 67; this 
would imply a level of chromospheric activity in the mid- 
dle between the Hyades and the old stars plateau, r ight 
in the Vaughan-Preston gap. lCarraro &: Chiosil l)l994|) . so 
far used as age reference, evaluate the age of NGC 752 
comparable, but slightly younger than IC 4651 (1.5 vs. 
1.6 Gyr). 

Before going into speculations about the difference 
between our and Barry et al.'s results, the NGC 752 
chromospheric-activity level should be proven by means 
of high-resolution spectroscopy in solar analogs. 

4.4. The impact of our results on the determination of 
ages and on the age of galactic population 

According to the picture we have drawn in the last Section, 
stars can be divided in three groups on the basis of their 
chromospheric-activity levels: 

- the youngest stars, with F^ > 1.2 ■ 10^ erg cm~^ sec~^ 

- the stars with F^ between ~ 0.7 and ~1.2 erg cm~^ 
sec~^ 

- the oldest stars, with F^ < 0.7 • 10^ erg cm~^ sec~^ 

Among our target stars, those in Praesepe and Hyades 
belong the first group, all the others to the third one. 
Let us call the upper age limit of the first ensemble <i 
and the lower age limit of second ^2- We don't know the 
exact values of ii and ^2- We only know that 0.6 Gyr 
< ii < t2 < 1.6 Gyr. 

In Figure [TUl we plot the histogram for 44 solar-type 
stars (spectral type bet ween GO and G 5) in the solar 
neighborhood taken from IPasauinil l|l992l) . The numbers 
of stars belonging to each group are respectively: 8, 2 and 
34. The result confirms what expected: most of the stars 
have an age older than ^2, much less are younger than 
ti, only ~ 5% have ages between ti and ^2, i-e. those in 
the Vaughan-Preston gap. In addition, the standard de- 
viation of the Fj- values of the 34 less active stars in the 



sample of lPasouinil ( 19921 1 is very similar to that of the 14 
targets belonging to the three intermediate age and old 

erg cm" 



clusters IC 4651, NGC 3680 and M 67: 0.12 10^ — --"^ 



see" . Unfortunately our results indicate that an accurate 
age determination based on chromospheric activity may 
not be obtained for stars older than ^ 1.5 Gyr. Once ti 
and ^2 are determined precisely, much could be said about 
the recent Star Formation History in the solar neighbor- 
hood, provided that initial conditions do not dominate the 
fast decline of chromospheric activity occurring between 
ti and ^2- 
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5. Conclusions 

We have studied the evolution of chromospheric activity 
and rotation in solar-type stars by analysing VLT and 
KECK high-resolution spectra of 35 stars belonging to 5 
open clusters and published solar data. As a measure of 
chromospheric activity we have used F^, i.e. the equiv- 
alent width of the Ca II K line emission core corrected 
for photospheric contribution and transformed into flux 
at the stellar surface. Our data set can be divided in three 
age bins: 6 Praesepe stars and 15 Hyades stars at 0.6 Gyr; 
5 IC 4651 stars and 2 stars in NGC 3680, at about 1.7 
Gyr; and finally 7 M 67 stars and the Sun, at about 4.5 
Gyr. We undertook the present research mainly to ascer- 
tain a decay law which allows measuring chromospheric 
ages in field stars. We have found that stars within each 
age bin show indeed the same mean level of chromospheric 
activity, confirming that it is possible to search for a chro- 
mospheric activity-age relationship. Stars in the youngest 
bin have chromospheric fluxes ranging between 1.2 and 2.8 
10^ erg cm~^ sec~^. All the other stars have values rang- 
ing between 0.3 and 0.7 10® erg cm~^ sec^^. The most 
important result is that there is no difference in the dis- 
tribution of chromospheric-activity levels between the in- 
termediate and the oldest age bin, even though difference 
of ages of the two bins is as large as ~ 3 Gyr. We obtain 
a very similar conclusion about the angular momentum 
evolution, with the only difference that the gap between 
the youngest bin and the others is less pronounced. A fast 
decay of chromospheric activity occurs roughly between 
0.6 and 1.5 Gyr of Main Sequence lifetime, after which a 
kind of plateau appears. 

No power law so far proposed can fit our data. The best 
match curve we computed, namely F^ = 0.5 • f^^'^ + 0.4, 
where F'j^ is in unit of 10®-erg-cm~^-sec^^, and the age is 
expressed in Gyr, can reproduce this trend. 

Coupling this trend with the observational fact that in 
each of the clusters observed the star to star activity level 
vary almost by a factor 2, our results indicate that the 
determination of chromospheric ages for field stars older 
than ^ 1.5 Gyr is extremely challenging, or even poten- 
tially misleading in particular if only one or a few activity 
measurements are available. 
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Table 1. I n this t able a ll of the meas urements reRarding our s ample stars are s hown. The na mes in Colu mn 1 are 
taken from lEggenllll969l) for NGC 3680 . lAnthonv Twarog et alJ lll98J^_(AMC ) and lEggenI lll97ll) for IC 4651. ISander^ 



l)l977|) for M 67, iKlein WassinkI l)l927J) for Praese pe and Ivan BuereiJ lll952l) for the Hyades stars. The photometry is 
taken froni the reference pap ers for IC 4651. from [ Nordstrom et al.l l|1997r for NGC 3680, from Hipparcos Catalogue 
|IeSAII1997I) for Hva des. from l.Tones fc Stauffeil lll99ll) for Praese pe. with the exception of KW 368 which is taken from 
I.Tones fc CudworthI (11,983), and from iMontgomerv et al.l l| 19931) for M 67. Regarding chromospheric-activity indices 
measurements, there are the 1-A indices derived from spectra either non corrected (Columns 4) or corrected (Columns 
5) for the absorption of the IS medium. To the corrected index we have subtracted the photospheric contribution, 
which is 0.049 A, and transformed the result into flux. The final data are in Column 6. v sini measurements in Column 
8 are taken from IPaulson et al.l 1(20031) for Hyades, for all the other clusters are based on the FWHMs in Column 7 
via the calibrations in Section 13.21 
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Fig. 9. The evolution of chromospheric activity (measured by F'j^ corrected for IS absorption, upper panel) and 
rotation (lower panel) with age. Here the vertical bars represent the whole spread of the ordinate data within the 
cluster. The points are the mean values. For both quantities the surprising result is that the intermediate age clusters 
show the same level of chromospheric activity and mean rotational velocity of M 67 and the Sun. Since, as remarked, 
the bars are not error bars but peak to peak spreads, their position attests that the activity distribution of the 7 
intermediate age stars and that of the 7 M 67 stars are virtually equal. For both diagrams we trie d to fit three power 
law s, as indicated in t he figure. -2/3 and -1/2 are the exponent values suggested respectively bv iBarrv et al.l (|l983) 
and lSkumanichl l|1972|) . As can be seen, they do not fit our data. By minimizing the x^ we find, for the free parameters: 
ai = -0.38, 6i = 1.79,02 = -0.79,62 = 2.24, 03 = 0.47,^3 = 0.43, exp = 2.79, ci = 0.71, di = 4.20, ca = -0.05,^2 = 
5.00, C3 = 2.23, dy, = 2.35, expl = 1.47. 
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Fig. 10. Histogram of the chromospheric-activity levels in the sample of stars from lPasauinil l|1992|) . 



